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Abstract 

We present some exact solutions of the Einstein equations with an anisotropic fluid 
exploiting the Chaplygin equation of state. The solutions describe spacetimes with two 
identical T regions and an intermediate static spherically symmetric R region containing 
a wormhole. The metric in the T region represents an anisotropic Kantowski-Sachs cos- 
mological model. Its evolution starts from a horizon and develops according to different 
scenarios including eternal expansion, contraction and also a finite universe lifetime. 
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1 Introduction 



Observational data of the last decade, such as Type la Supernovae (SNIa) [1], Cosmic Microwave 
Background (CMB) [2] and Large Scale Structure [3], are evidence of accelerating expansion of 
the Universe. This behavior is caused by the dark energy featured by a negative ratio of the 
pressure to the energy density. A dark energy nature is not clear nowadays and a number of 
competing models describe the dark energy as an exotic matter or as a result of non-minimal 
gravitational interaction with ordinary matter. 

The simplest dark energy model, the cosmological constant, is indeed the vacuum energy 
with the equation of state p = —p. A number of models, such as quintessence [4], k-essence [5], 
phantom [6] and etc., are based on scalar field theories. An alternative way is to be thought 
a geometric origin of the dark energy. It motivates to modify the four-dimensional gravity 
theory [7] or to introduce additional dimensions as in braneworld models [8]. The Chaplygin 
gas model, also denoted as quartessence, exploits a negative pressure fluid, which is inversely 
proportional to the energy density p = —a/p [9]. It was modified for the models unifying the 
dark energy and dark matter [10], including generalized Chaplygin gas [11, 12] and superfluid 
Chaplygin gas [13]. For a more detail review of dark energy models and references see [14]. 

In the cosmological context the Chaplygin gas originally appeared as a hydrodynamical 
representation of the generalized Born-Infeld Lagrangian 



describing a (3+l)-dimensional brane universe with the scalar field 9 in a (4+l)-dimensional 
bulk [11]. It was also shown that the Chaplygin gas model can be derived from the Lagrangian 



for a complex scalar field in the WKB-approximation [13]. 

The Chaplygin equation of state can violate the strong energy condition p + 3p > 0, which 
is necessary for an accelerating cosmological process as well as it can violate the null energy 
condition p + p > providing for the existence of traversable wormholes. Wormholes with 
the Chaplygin gas are originally studied by Lobo [15]. Owing to the peculiar equation of 
state the obtained solutions possess some features. Primarily, if the solution exists in the 
whole spacetime then it cannot be asymptotically flat. To ensure asymptotic flatness for the 
Chaplygin wormhole the interior solution in [15] is matched with the exterior Schwarzschild 
solution using the Darmois-Israel formalism of a dynamic thin shell [16]. This method is also 
used when the wormhole solution exists only within a spatial sphere and cannot be directly 
continued outside because it is singular on the sphere. 

An alternative approach to the Chaplygin wormholes is to consider horizons instead of 
singularities. A cosmic expansion beyond an event horizon can be considered as a model of an 
origin of our Universe [17, 18, 19, 20]. The idea that the universe can emerge in the interior of a 
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black hole was proposed by Pathria [17]. The primary attention to the universes inside the black 
holes is concentrated on the inquiry of geometrical properties and cosmological details rather 
than on exact solutions of the Einstein equations [18]. Bronnikov et al. [19] used the exact 
solutions for phantom scalar fields to study and classify regular black holes where singularity 
is replaced by a cosmological expansion and which therefore were called black universes. Such 
objects can be supported by phantom matters and could arise from collapse in another ambient 
universe. In the work [20] the Lemaitre type cosmology for an anisotropic perfect fluid with 
the vacuum equation of state was considered for different topologies of the spacetime. The 
conditions for the cosmological evolution starting from a horizon for different kinds of matter 
were studied in [21, 22]. 

We apply the similar structure to Chaplygin wormholes. It is free of singularities and the 
wormhole is surrounded by horizons separating the spacetime into the static and cosmological 
regions. In this work we deal with the simple horizons for which we determine the conditions 
of regularity for the metric and the Chaplygin matter. Some exact solutions corresponding to 
different cosmological scenarios are developed using these conditions. 

The paper is organized as follows. In the next section we consider the general properties 
of spacetimes with a spherically symmetric wormhole surrounded by two horizons and external 
universes which are supported by matter with the Chaplygin equation of state. Some particular 
cosmological scenarios are presented in Sec. 3. Concluding remarks are summarized in Sec. 4. 



2 Wormhole between two horizons 

We consider a static spherically symmetric spacetime with the metric 

ds 2 = Air) dt 2 ^Tr" - r H d # 2 + sin 2 9 dip 2 ) (3) 

1 — b{r)/r 

and the stress-energy tensor of an anisotropic perfect fluid in which the energy density and the 
radial pressure are related by Chaplygin's equation of state 

Pr = , (4) 

P 

where a is a positive constant. In this case the Einstein equations = StcT^ (we set 
c — G — 1) reduce to 

b' = 8nr 2 p, (5) 
A' b + 87rr 3 p r 

~A = r 2 (l -b/r) 1 ( ' 

2 A' 

Pr = -(Pt-Pr) ~ (P+Pr)~£, (7) 

where the prime denotes the derivative with respect to r and pt(r) is the transversal pressure. 
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The metric (3) describes a static spherically symmetric wormhole if the radial coordinate r 
increases from a minimum value, corresponding to the wormhole throat, to infinity. Without 
loss in generality we take r m j n = 1. In this approach r is associated with two maps covering 
the spacetime on either side of the throat. The functions A(r) and b(r) are regular and defined 
through all the range 1 < r < oo and satisfy the flare-out conditions 

6(1) = 1, fe'(l) < 1 (8) 
0<A(l)<oo, A'(l)<oo (9) 

to provide a minimal area for coordinate spheres at r = 1 and the metric regularity on the 
throat. 

For traversable wormholes it is also necessary that A(r) should be positive and b(r) < r 
elsewhere to avoid horizons and singularities. These requirements are ruled out in the present 
paper in favor of the horizon located at the radius r = r > 1. 

The horizon conditions correspond to simultaneous solutions of the equations b(r ) = r 
and A(ro) = as it is depicted in Fig. 1. We deal with only simple horizons. It means that the 
metric component g rr has a simple pole at the point r , so that 

b(r) = r + b'(r )(r - r ) + o(r - r ) (10) 

and the function A(r) has a simple zero 

A(r) oc (r — r ) + o(r — r ). (11) 

These horizon conditions and compatibility of the Einstein equations near the horizon lead to 
the following relation for the Chaplygin gas 

b'(r ) = 8nr 2 ^. (12) 

In this case r corresponds to the simple horizon separating a static R region at r < r from a 
non-static T region at r > r . Indeed, for the metric (3) the square of the normal to a surface 
r=const is A = bjr — 1 and A < when 1 < r < r®. In this region the spacetime is static 
(R region) and r is a spatial coordinate. At r > r the normal is a spacelike vector and A > 0. 
This is a non-static T region where r becomes a temporal coordinate. 

This class of metrics is regular in any finite region for r > 1. It is verified using the 
Kretschmann scalar K = R^sR 01 ^ 6 , where R a /3j8 is the Riemann curvature tensor. To avoid 
too cumbersome representation for the Kretschmann scalar we write it in the form 

_ P(r,A,A>,A»,b,l/) 

where the numerator is a polynomial function depending on the indicated arguments and P oc 
(r — r ) 4 + o([r — r ] 4 ) near the horizon r = r . Taking into account the expansion (11) it 
means that the Kretschmann scalar is finite on the horizon. 
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Figure 1: The functions Air) and b(r) for 
the geometry with a wormhole at r = 1 
and simple horizons at r = r . 



Figure 2: The global structure diagram of the 
model. The thick lines correspond to r = oo 
and the thin lines correspond to the horizons at 
r = r . 



The metric (3) in the T region can be rewritten as 

d, 2 = ^- I -i(,)dp 2 -, 2 dfi 2 , (14) 

where A(r) = —A(r) is the positive function when r > Tq. We introduce the new coordinates, 
77 instead r and p instead t to mark the interchange of roles between the spatial and temporal 
coordinates. The integral 

r = ±[ - dV (15) 

J VKv)/v-i 

defines the physical time r for an observer in the T region. The metric (14) describes a 
Kantowski-Sachs (KS) anisotropic cosmological model [23] with two scale factors for the spher- 
ical and longitudinal directions which depend on the proper time r according to a f (r) = r/(r) 
and a r (r) = A 1 / 2 (77(7-)) correspondingly. The global structure of the spacetime is shown in 
Fig. 2. The metric (14) describes both the contracting region T_ corresponding to r < and 
the expanding region T + corresponding to r > 0. 

At r) = To there is a coordinate singularity which the observer interprets as the beginning of 
the universe. In [20] it is called a null bang sinse a horizon is a null surface where the spatial 
volume vanishes. It is easy to see from Eq. (15) that if the function b — > 00 faster than if 
for 77—7-00, then the universe evolves during a finite cosmological time and otherwise it lives 
endlessly. 



It is evident from the Einstein equations that the following conditions are fulfilled on the 
horizon 

p(r ) = -p r (r ) = y/a, (16) 

that corresponds to the requirement p r + p = on the horizon [21] which is independent of 
the matter kind. The condition (16) provides continuity of the energy density and the radial 
pressure in going from R to T region. In fact, p and p r interchange their positions in the stress- 
energy tensor when crossing the horizon. The equation of state (4) is precisely the same in 
both regions. Under the conditions (11) and (12) the transversal pressure is also a continuous 
function and takes the value Ptij'o) = b"(r )/8irro — on the horizon. 

Below some exact solutions of the equations (4)-(7) with the conditions (8)— (12) corre- 
sponding to different cosmological scenarios are studied. 



3 Different types of Kantowski- Sachs cosmologies 

The equation of state (4) relates the energy density with the radial pressure only. It allows 
one to solve the problem predetermining the metric function b(r) satisfying the flare-out and 
horizon conditions. The other functions are found from Eqs. (5)-(7) and produce different 
cosmological scenarios in the T regions whereas the R regions are similar in all solutions. 

3.1 Eternally expanding universes 

One of the simplest functions satisfying the conditions (8) and (10) is the parabola 

b( r ) = d(r - l)(r - r ) + r. (17) 
The function A(r) is expressed from Eq. (6) with the conditions (9) and (11) as 

A(r) = A {r °- r){r + qi)q2 e^\ (18) 



(19) 



where 

r d + d - 1 2q$ 1 - d + r d ( 1 

qi= 2d ' q2 = ^f l + d-r^ * = 2 *{** + d 
_ (r d-d+l) 2 ,_ r 2 y^T8-r 4 -2 
q4 ~ Kd 2 ' 2(r -l) 
are positive constants and r is connected with the parameter a from the equation of state (4) 
by the relation 

(1 - d - r Q d? ( . 

In the R region the metric (3) can be transformed to the form 

ds 2 = B(x) dt 2 - dx 2 -r 2 (x)(d6 2 + sin 2 9 dip 2 ), (21) 
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Figure 3: The metric functions depending on the 
Gaussian coordinates in the R region. This fig- 
ure is plotted for b(r) given by Eq. (17) when 
ro — 1.5 and other parameters are calculated ac- 
cording to (19), however the functions r(x) and 
B(x) have a similar form for the other solutions. 

where 

r(x) = 1 + (r - I) sin 2 e (yd x/ , (22) 

and e(z) is the inverse function to the elliptic integral of the second kind -E^l — T'o)- The 
metric function B(x) = A(r(x)), where r(x) from Eq. (22) is substituted into the expression 
(18). Behaviors of r(x) and B(x) are depicted in Fig. 3. 

Contrary to the representation (3), the metric (21) describes the whole R region with 
the Gaussian radial coordinate x lying in the range \x\ < irE (n/2\l — ro)/Vd- When x — > 
±nE(n/2\l - r )/Vd the function B (x) becomes zero providing the evidence for horizons 
which are spaced at a finite distance apart for a static observer. The function r(x) has a mini- 
mum at x = implying a wormhole throat, so the R region is a spherically symmetric one with 
the wormhole surrounded by two horizons. 

For observers in the T region the universe evolution starts from the horizon at some instant 
in the past and lasts for infinite time. All this time the universe expands in all directions as 
it demonstrates in Fig. 4. In the far future the scale factor along the coordinate spheres is 
quadratic in the proper time r while the scale factor in the radial direction grows exponentially 

a r oc r q2 exp g 4 T 4 /2, a t oc r 2 . (23) 
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Figure 4: The increase of the scale 
factors for the solution specified by 
Eqs. (17)-(19) for r = 1.5. 



Figure 5: The dimensionless energy density, ra- 
dial and transversal pressures as functions of the 
radial coordinate for the eternally expanding uni- 
verse. 



The behaviors of the energy density and the pressures are shown in Fig. 5. The energy 
density is positive and the radial pressure is negative in the whole spacetime whereas the 
transversal pressure reverses its sign. In the T region all these quantities are comparable close 
to the null bang. At the late stage the radial pressure decreases in its absolute value as p r oc t~ 2 
and the main contribution is given by the energy density growing quadratically in the proper 
time, p oc r 2 , and the negative transversal pressure which increases in its absolute value as 
p t oc r 6 . 



3.2 Universes with finite lifetime 

The alternative scenario is arrived from the equations (4)-(7) if the function 

b(r) = dr A + l-d. (24) 

In this case the function 

Mr) = 'T^ -P (-^1 arctan ( *L±2±1) ) , ( 25 ) 



r[r 2 + (r + l)r + qjq 2 + r ] qi { q 2 \ g 2 

where 

Qi = ^^, Q2 = ^ + 2^ + 3, r = l + V2. (26) 

This solution implies a = (3 — 2\/2)/647r 2 in the equation of state (4). 

In the R region the spacetime has the same structure as before. The metric can be trans- 
formed to the form (21) with the only difference that the metric functions B(x) and r(x) cannot 
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Figure 6: The behaviors of the scale fac- 
tors for the solution when the function 
b(r) and A(r) are given by Eqs. (24) and 
(25). The proper lifetime of the universe 
is finite. 



r 

Figure 7: The energy density, radial and transver- 
sal pressures as a functions of the radial coordi- 
nate in the universe with the finite proper lifetime. 



be obtained analytically and all calculations are numerical. Their qualitative behaviors are the 
same as in Fig. 3 

The expansion in the T region starts with a null bang and the ensuing cosmological evolution 
develops during finite physical lifetime from the point of view of the KS observer. Fig 6 
demonstrates that throughout this period the spherical scale factor infinitely increases and 
the longitudinal scale factor peaks at some instant and falls to zero, although both radial 
and transversal pressures are negative and the energy density disappears towards the end of 
the evolution as it is shown in Fig. 7. This situation arises in the KS model from the fact 
that a transversal acceleration is determined by the radial pressure and the transversal Habble 
parameter H t = a t /a t , where the dot denotes the derivative with respect to the proper time 
r, while a radial acceleration also depends on the transversal pressure and the energy density 
according to 

\ + Hi + 2^ = -87rj9 r , (27) 
J_ + E 2 t - ^ = 4tt (p - Pr + 2jh) . (28) 

If the energy density is neglected and the rate of growth of H t is small so that \p r \ > 2\pt\ + 
Hi /4vr, then one observes a decelerated expansion or contraction in the longitudinal direction. 
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3.3 Eternally contracting universes 

An intermediate scenario is obtained for the metric functions 



where 



b( r ) = dr(r - l)(r - r ) + r, (29) 

Air) = A — F Tororexp < —2 arctan > , 30 

^ \^[3^-2 qi r+r q 3 ] q2q5 P \ q 6 \ q & ) J' K ' 



qi=r Q + 1, q 2 = d{2 + 5r + 2r 2 ) - 2d 2 r (l + r + r 2 ) - 1, 
g 3 = 1 + ^r, ?4 = d(2 + 5r + 2r 2 ) - d 2 r (2 - r + 2r 2 ) - 4, 

n. - 1 - ^ + ^ _,/CTZ^ ( 31 ) 

6rf 2 r4(l + rf-rfr )' m y_d 



(1 - rfr + rfr 2 ) 2 r 0V /r^ + 4r 2 + 4r -^-2 

647r 2 r 2 ' 2(r 2 -l) 

are positive constants except q 2 , which can be both positive and negative 1 , and the horizon 
radius lies in the range 1 < r < r max 3.04. If r > r max the function A(r) resulting from 
Eq. (6) gets singularities in the static region. 

Note that the transformation to the metric (21) in the R region gives the quite simple form 
for spherical radii 

r(x) = i (Vo + 1) - (r - 1) cos Vdx^ (32) 

with the Gaussian coordinate x lying in the range [—n/Vd,7r/Vd\. 

For the KS observer the universe starts to evolve from the horizon and gradually dies into 
infinite future because the expansion in the radial direction gives place to the contraction. The 
scale factors for the longitudinal and lateral directions as functions of the proper time are shown 
in Fig. 8. The radial pressure is everywhere negative (see Fig. 9), so that it guarantees the 
accelerated expansion along the coordinate spheres, and the inequality p—p r + 2p t — H 2 /An > 
provides the decelerated growth and the subsequent decay of the radial scale factor according 
to Eq. (28). At the late stage the scale factors vary in the proper time r by the exponential 
laws 

a r oc exp <^ -Vd I q 2 q 5 - ttt- r , a t oc exp Vdr. (33) 

The asymptotic expression (33) resembles in some way the de Sitter behavior bringing the 
pressures and the energy density to the constant values 

8™ 3d d rj .„ \ / -i ^ m /„„x 

P ~*~M' Pr_> ~87' Pt ^ ~^^ + ^~ qi ~ 2q2qs ^ 1 ~ qi ~ 2q2q ^- ( 34 ) 

Such a tendency takes place for any function 6(77) increasing as rj 3 when r\ — > 00. It can provide 
the basis for a cosmological scenario with an asymptotic de Sitter phase at late times. 



1 Nevertheless the value q$ + 2q 2 q5 — 1 > for all values of r lying in the range 1 < r < r max . It means that 
A(r) vanishes when r — > 00. 
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Figure 8: The behaviors of the scale fac- 
tors in the solution with the metric func- 
tion (29) and (30) where the quantities 
(31) are calculated for r = 1.5. The 
universe expands in the spherical direc- 
tion and contracts in the radial direction 
during an infinite time. 

4 Conclusion 

We have considered some exact solutions of Einstein's equations with anisotropic fluid describ- 
ing a spacetime with two identical cosmological T regions and an intermediate static spherically 
symmetric R region. The latter contains a wormhole allowing one to pass between two horizons. 
The T regions are described by anisotropic KS models with different kinds of cosmological evo- 
lution. We have presented the solutions describing eternally expanding universes, disappearing 
ones and universes with a limited lifetime. 

These solutions rest upon the flare-out conditions (8), (9) and the horizon regularity con- 
ditions for matter with the Chaplygin equation of state (10)— (12), (16). Similar regularity 
conditions were derived before for the null bang cosmologies with different types of matters 
such as the neutral and charged matter with a barotropic equation of state p r ~ p m , m > 1 
[22] and matter with a non-vacuum behavior (i.e. the ratio of the radial pressure to the energy 
density w = p r /p ^ —1 everywhere) [21]. However they are incompatible with the model under 
consideration since the regularity for the Chaplygin gas is realized only through the "vacuum" 
value w = — 1 on the horizon. 

The presented solutions violate everywhere the dominated energy condition, p > 0, p > 
\Pr\, P > \Pt\ while the strong energy condition proves to be realized in some parts of the 



Figure 9: The energy density, radial and transver- 
sal pressures as functions of the radial coordinate 
in the gradually dying universe. 
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spacetime as it takes place in the last scenario. One sees from Fig. 9 that the relation p + p r + 
2p t > holds in the layer near the horizon. 

Note that the suggested spacetime conformation does not support an isotropic solutions 
in both R and T regions. Indeed, the isotropy of a spacetime implies p r = p t and Eq. (7) is 
readily integrated. The resulting function A(r) vanishes concurrently with the energy density 
and hence with the derivative b'(r). Fig. 1 pictorially shows that such a solution is incompatible 
with the considered structure. 

Besides, there is no an isotropization in the solutions obtained, although it has considerable 
utility for the early Universe model. Nevertheless, isotropization can be most likely realized 
in the studied model. The last scenario has the asymptotically constant pressures and energy 
density, and solutions of this kind may quite probably demonstrate the de Sitter behavior at 
late times for the special selection of the metric functions. This problem is the subject of our 
current studies and the results will be presented in our future papers. 
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